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Abstract. We establish contact between the delocalization properties of pure 
quantum states, as quantified by their number of principal components, and the 
average generalized entanglement properties, as quantified by purity measures 
relative to different observable sets. We find that correlations between products 
of state vector components with respect to Hamming distance play an important 
role in the structure of subsystem-based purity measures. In particular, we derive 
general conditions under which the amount of global multipartite entanglement 
relates to the inverse participation ratio averaged over a maximal set of mutually 
unbiased product bases. Furthermore, we provide a method for computing 
the expected amount of generalized entanglement with respect to an arbitrary 
observable set for random pure states. Specific examples and an explicit 
application to a disordered quantum spin chain are discussed. 
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1. Introduction 

Developing a deeper qualitative and quantitative understanding of "complex" 
quantum systems is a broad challenge whose implications range from condensed 
matter physics to fundamental quantum theory and quantum information science 
(QIS). In a loose sense, complexity may be intuitively associated with the lack of 
a "simple" description of physical properties in situations where such a description 
should in principle follow from a small set of known, basic rules .1 . In quantum 
systems, complex quantum features so defined may appear at both the kinematical and 
dynamical level via three main pathways: large state-space size; interaction between 
constituent subsystems; absence of dynamical regularities and non-integrability. Taken 
together, these factors may be ultimately held responsible for non-scalable (typically, 
exponentially inefRcient) parameterizations of system properties; the emergence of 
many-body phenomena and quantum irreversibility (via interactions both within the 
system of interest and between the system and its environment); and the possibility 
of dynamical instability and quantum chaos. 
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Physically, the occurrence of genuinely quantum correlations - entanglement - 
lies at the heart of all the above phenomena. From a QIS perspective, entanglement 
is also intimately tied to the inherent complexity that distinguishes quantum from 
purely classical information processing - being, in particular, the defining resource 
for quantum communication [2] as well as a necessary prerequisite for quantum 
computational speed-up [3]. In recent years, this has naturally motivated a host 
of investigations aimed at characterizing the nature and role of entanglement in 
complex quantum systems. While a complete understanding is far from being reached, 
important progress is being made toward elucidating the behavior of entanglement 
across a quantum phase transition [4l [5] and across a transition from integrability 
to quantum chaos [3 [7l [8], as well as in better accounting for entanglement in 
computational schemes for interacting quantum systems, like renormalization group 
methods 9 . 

In this context, the notion of generalized entanglement (GE) [10] has recently 
emerged as a unifying framework for describing entanglement in arbitrary physical 
and QIS settings ~ capable, in particular, to recover conventional (subsystem- 
based) entanglement in well-defined conditions and to directly incorporate physical 
constraints such as quantum indistinguishability. Beside providing new Lie-algebraic 
measures for diagnosing broken-symmetry quantum phase transitions in a variety of 
models [S], GE has contributed so far to the understanding of standard multipartite 
spin correlations in disordered lattice systems [IT], provided a natural testbed for 
investigating entanglement generation in chaotic quantum maps |12j . as well as shed 
light on conditions for quantum-computational speed-up in a wide class of Lie- 
algebraic models [13j. 

Here, we continue to explore applications of the GE framework to complex 
quantum systems by focusing, in particular, on highlighting the relationship between 
GE measures and state delocalization properties^ as quantified by standard indicators 
like the number of principal components with respect to an appropriate basis. In 
line with the growing body of work at the interface between QIS, condensed-matter, 
and quantum statistical physics [14j . our main motivation is to qualitatively and 
quantitatively characterize points of contact between notions originally developed 
in different contexts - in the hope that this may lead to useful cross-implications. 
Following a review of the essential GE background in Section 2, the relationship 
between global multipartite entanglement and delocalization is addressed in Section 3, 
by explicitly uncovering the role of Hamming distance in the structure of subsystem- 
based GE. In Section 4, a general method for estimating the expected GE of random 
pure states is presented. In Section 5, a concrete application to a disordered Heisenberg 
spin- 1/2 chain is discussed. Final remarks conclude in Section 6. 

2. Generalized entanglement and purity measures 

The basic idea of the GE approach is that the entanglement properties of quantum 
states are determined by the expectation values of a distinguished subspace of 
observables rather than a preferred decomposition of the system into subsystems 
[TU] . This allows the notion of GE to be meaningful in physical settings which 
the conventional notion is too narrow to embrace: in particular, GE is directly 
applicable to systems subject to limitations in the available control interactions and 
measurements, and described by arbitrary operator languages (spin, fermion, etc) 
as for instance many-body quantum systems. In addition, because the GE notion 
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rests only on convexity properties of spaces of quantum states and observables, GE is 
mathematically suited for entanglement formulations in abstract operational theories. 
We refer the reader to (TUl [HI HH HZl HH] for a more expanded discussion. 

The starting point for defining GE is to realize that a pure state of a composite 
quantum system is entangled (in the usual sense) iff at least one of the reduced 
subsystem states is mixed. Let the system of interest be described by a pure state, 
\tp) & H , with H, dim(H) = N, and p = being the associated Hilbert space 

and density operator, respectively. Let in addition the distinguished observable set 
consists of the Hermitian operators in a linear subspace h C B{H) of the full operator 
space on 7i, with h closed under Hermitian conjugation. The key step is to replace 
the notion of reduced state as obtained via a partial trace in the usual tensor-product 
sense by a notion of reduced state as resulting from the restriction to h of the positive 
linear functional to corresponding to p via the trace map [lOl [TBI [E] • Such a reduction 
may be specified in terms of the (unique) projection map Vh with respect to the trace 
inner product, p i— > Vh{p)- Accordingly, is defined to be generalized unentangled 
relative to h iff its reduced state Ph{p) is pure - that is, extremal in the space of 
reduced states [19] . 

While no unique measure suffice to quantify the amount of /i-GE present in jV'}, 
the simplest possibility is to evaluate the square length of Vhip)- Let be a 

basis of Hermitian traceless operators for h, orthogonal in the trace inner product, 
tr(feifej) = NSij. The purity of \ip) relative to h (/i-purity) is given by 

PiM) = n^J2^r{phf = A^,^ (1) 

i i 

where the normalization constant Hh depends in general on h and N, and ensures that 
the maximum value of Ph is 1. Thus, a state [ip) with maximal purity, Phdijj)) = 1, 
is unentangled with respect to h, hence it has extremal length [17 . In the physically 
relevant case where h forms a (irreducibly represented) Lie algebra, maximal /i-purity 
is both a necessary and sufficient condition for a pure state to be /i-unentangled, Ph 
is invariant under unitary transformations generated by arbitrary elements of h, and 
GE/j = 1 — is an entanglement monotone [10 . Ph may be extended to a measure 
for mixed-state GE via a standard convex-roof construction [10] . 

The following specializations and applications of the above GE definition may 
serve to clarify the relationship between GE and standard entanglement, and will be 
especially relevant for the present discussion: 

(1) Absolute purity. By definition, p is pure iff it is a one-dimensional projector, 
hence iff tr(p^) = 1. By identifying h with the (real) Lie algebra of all traceless 
observables onTi., h = su{N), Eq. ^ gives Kaii = 1/{N ~ 1) and 

Pauim = P.uiN){m - J^{^<P') - (2) 

consistently normalized so that purity is 1 for pure states and for totally mixed ones. 

(2) Bipartite systems and linear entropy. For a system consisting of two 
subsystems A and B,Ti. — 'Ha®'Hb, dim(7-^yi) = dA, dim(7iB) — ds, the information 
accessible through measurements on A or _B alone is contained in the reduced density 
operators pA = trB(p), pB ~ tr^(/9). As mentioned, a pure state {i/j) G H is 
unentangled, \ip) = \4>a)'Si\4>b) , iff both pA and pB, are pure. A bipartite entanglement 
measure known as the linear entropy E (of either subsystem) may be constructed 
as EaUiP)) = 1 — tr(/9^). In the GE approach, standard bipartite entanglement is 
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recovered by choosing the set of all uni-local observables acting on A (or B) alone, 
e.g. h^JiA^ Bu(dA) © Ib (equivalently, e''' = SU(dA) ®^b)- Eq. (P) then yields 

1 - PnAm = 1 - rff^(tr(^A) - ^) = rffri^^(l^))' (3) 

that is, GEifiA directly proportional to the linear subsystem entropy. 

(3) Multipartite systems and average subsystem purity. The above example 
generalizes to a multipartite system consisting of n subsystems of dimension d. That is, 
conventional (subsystem-based) entanglement is recovered by selecting the algebra of 
all uni-local observables acting on individual subsystems as distinguished observables, 
h = hioc = ®ihi = sui(d) ® . . . © su„(d). Eq. ([T]) then gives 

that is, the purity with respect to arbitrary local observables is equal to the average 
(normalized) subsystem purity, as intuitively expected ,16j. Pioc{\i^)) attains its 
maximum 1 only for completely separable states, 1-0) = and is equal to iff 

each reduced density matrix is totally mixed (hence no information is available through 
local operations). The entanglement measure GE/oc = 1 — Pioc is thus proportional 
to the average linear entropy over all bi-partitions of the system into blocks of 1 and 
(n — 1) subsystems. For qubit systems [d = 2), such a measure has been shown in 
[5] and |20j to coincide with global multipartite entanglement Q as introduced by 
Meyer- Wallach [gT], Q(|^)) = 1 - F;oc(|V'))- 

(4) Expected h-purity of a set of states. For a fixed observable set, the expected 
/i-purity of a pure state taken with respect to a certain probability distribution ^ 
quantifies GE properties of a typical state in the ensemble, 

P„=E(«){P;,(|7/^))}. (5) 

An important instance arises for uniformly sampled random pure states, in which case 
^ coincides with the unitarily invariant Haar measure on SU(7V) |22j . 



3. Delocalization and local purity 

Given an orthonormal basis {\k)} in H, a well-established measure of state 
delocalization in quantum statistical physics and quantum chaos is the number of 
principal components (NPC), 

NPC(|V')) = (El(^l'/^)lT'^(^l"'=lT'' El«'=l' = l' (6) 

k k k 

quantifying the number of basis states on which lip) has a significant amplitude G C. 
NPC so defined ranges from a minimum value of 1, meaning that \^p) coincides with a 
single basis element, to a maximum of A^, corresponding to a maximally delocalized 
state with equal probabilities ja^p = 1/N. If NPC > 1, measurements in the 
corresponding basis will result in a probability distribution over possible outcomes. 
For instance, a crossover from localization to delocalization with respect to a large 
number of basis states occurs in the eigenvectors of the Anderson model during the 
insulator-to-metal transition, as well as in the eigenvectors of quantum spin lattices 
across a transition to quantum chaos [6J . NPC is equivalently referred to as participation 
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ratio (or participation number [7]). Accordingly, NPC^^ will be often denoted here as 
inverse participation ratio (IPR) [231 124] . 

From a conceptual standpoint, it is interesting to observe that NPC (and IPR) may 
be directly related to an appropriate /i-purity. Specifically, let h = hdiag denote the 
subspace of all (traceless) observables which are diagonal in the chosen orthonornial 
basis {\k)}. Then, 

N 1 1 

Ph,..A\^)) = iV-lNPC(l^)) " TV^- 

Such an observable space may be considered the commutant of a non-degenerate 
Hamiltonian and forms a trivial, abelian Lie algebra. As such hdiag does not identify a 
decomposition into quantum subsystems, and G^diag need not have any relationship 
to entanglement in the standard sense. Clearly, PhdiagiW) = 1 ifi^ NPC(|-;/;)) = 1. 
Thus, in a sense, Phdiagi\''P)) may be also thought as quantifying how non-classical 
is relative to the given basis. 

Our next objective is to investigate to what extent a relation similar to Eq. ([7]) 
may exist between standard entanglement (as quantified by Pioc) and NPC, as evaluated 
in each of a maximal set of mutually unbiased product bases [25[ . A product basis 
is one where each basis state is unentangled. Two bases are mutually unbiased if 
localization in one basis implies maximal delocalization in the other. In general, we 
shall find that Pioc is not solely a function of NPC in these bases, but it also depends 
on additional structure of the input state. 

Focus on a system consisting of n qubits (spin-1/2) first, N — 2". The bases 
{|/cq)}, consisting of the joint eigenstates of qubit observables {ctq''}, a — x,y,z, 
provide a natural maximal set of mutually unbiased product bases. Let {a^} denote 
the components of | V'} in the basis { Ifco.) } . The local purity Pioc may then be expressed 

as Pioc = Px+ Py+ Pz, where Pa(IV')) = ^ V)^ [H [S]. Recall that the 

Hamming distance between two binary bit strings of equal length measures the number 
of substitutions required to change one into the other. We have: 

Lemma 3.1. For every pure state j?/"} of n qubits, the following identity holds: 
P„(|^)) = l-l5]A.,Kp|a^"p, (8) 

where fkj is the Hamming distance between basis states \ka) and jja) that is, the 
number of instances where the eigenvalues of CTq'' differ on \ka) and \ja)- 
Proof. Note that we may express 

k k' 

where the (un)primed sum is over all such that the fc-th basis state has a (0)1 for 
the i-th qubit. Squaring both sides yields 

{^w^:^w = J2k\'+T,K'\' (9) 

k k' 

I o/'X^I a\2\ a\2 , \ ^ i a\2\ a \2 I a\2\ a \2\ 

+ ^[2^\"-k\ \aj\ + 2^ Iflfc'l Iflj-'l -Z^lflfcl \ak'\ )■ 

k<j k' <j' kk' 
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From the normalization of ji/i) one obtains 

5]|a^|'*^IPR„ = l-25]|a^na^"p. (10) 

k k<j 

Substituting ^ into ^ yields 



(^|<|V')^ = l-4V|a^na^,| 



2 



kk' 



Note that this sum is over all pairs such that \ka) and \k'^) differ on the ith qubit. 
Hence, the number of occurrences of a term |a"p|a"p in the sum yielding P^, which 
is over all qubits, is equal to the Hamming distance between \ka) and \ja), whereby 
the result. □ 

Note the structural similarity between the expressions for Pa and IPRa, Eqs. ([8]) 
and pO|) ~ the main difference being that in Pa the products |a^p|a"p are weighted 
by Hamming distance whereas in IPRq they are not. We have the following: 

Theorem 3.1. Assume that for each basis a = x,y,z, the values of the 
terms |afcP|a"P are independent on average upon Hamming distance (uncorrelation 
assumption). Then for every pure state of n qubits, 

where = 2" is the dimension of the Hilbert space. 



Proof. Let AJ — |aj?p|a"p denote the average over all pairs k,j, constrained 

to a specific Hamming-distance value fkj = f, and let A" = |a^P|a"|2 denote the 
unconstrained average over all k and j. Then the weighted average over all pairs k,j 
may be separated into the sum of averages over pairs corresponding to a given /, 

k,j f f f 

where is the number of pairs k,j with fixed Hamming distance /. Under the 
uncorrelation assumption, each of the ratios — 1, irrespective of /. Thus, by 
invoking the expression of Pa in Lemma 3.1, and by making the average over pairs 
defining A" explicit, 

4/,-- . „A / 2 



^"(i^» = i--(E-//-/)( 



f ' ' k<j 

4 



nN{N 



-T)(E"//)(i-™"(l^)) 



/ 

To evaluate J^f^ff^ ^^st note that for each state, |fc), there are (Jj) states labelled 



by j that are Hamming distance / from |fc). Thus, n/ = (y). Using = („"j), it 



N /n^ 

. LJ SI 

follows that E';=o/(;) = tE/=o (;) = t^- Hence, E/"// = t^- By summing 
over a, the result follows. □ 

Thus, Pioc depends in general on both the NPC in a set of three mutually unbiased 
product bases and on the average correlation of the products |afcP|ajP with respect to 
Hamming distance in each basis. 
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3.1. Conditions for single-basis delocalization 

For states obeying certain symmetries, the number of bases involved in the relationship 
between delocalization and Pioc may be reduced. 

A first physically relevant example is provided by states invariant under a non- 
standard (anti-unitary) time reversal symmetry T such that = I [23]. All states 
invariant under T may be expressed using only real components in an appropriate 
basis. For such states the expectation values involving the imaginary part of the 
operator space are zero. For instance, for states of qubit systems which are real 
in the standard (Ifcz)} basis, it follows that {ip\ay'' {iJj) — for all i. Hence, Py = 0, 
and NPCj^ does not enter the expression for Pioc. 

Notably, a further simplification occurs for the energy eigenstates of a large class 
of two-body spin Hamiltonians, which includes the Heisenberg, XXZ, XY, and Ising 
models - specifically, any Hamiltonian which may be written in the form 

for arbitrary coupling parameters Ja G M and on-site energy splittings Si G M. 
Any such Hamiltonian commutes with the collective Pauli operator (^"^^ cri*\ which 
describes a global Z2 symmetry. If H is non-degenerate^ then each eigenvector is 
invariant under this symmetry. It then follows from standard properties of Pauli 
operators (namely, that if [iTa,crf,] ^ and UalV') — -^alV'): then {ip\ab\fp) — 0) that 

(V^|a«|V') = (V|a«|V) = 0, Vz. (12) 

Under such conditions, Pioc = Pz, hence global entanglement properties, depend only 
on NPC2 and Hamming-correlations in the {Ifcz}} basis pS] . 

Note that if = , Vi , then (^"^j^ o'x'' and (^"^j^ cry*'' are also symmetries of H, 
resulting in (■0|cri*'' jV'} =0, Vi , thereby yielding Pioc = for each eigenvector of H. 

An important class of states obeying Eq. p2|) are the eigenstates of total z-angular 
momentum, Sz = ^^^"^ ■ For the Sz = subspace, for instance, the assumption of 
no correlation between Aj and fkj implies 

^iS^=0)fUM _ ^0 1 L 

iVo~lNPC,(|V)) iVo-1' 

where Nq = n!/[(n/2)!]2 is the dimension of the subspace. In the n-dimensional 
Sz = n — 2 subspace describing the single-excitation sector, fkj = 2 for all pairs of 
basis states, hence Pioc depends directly on NPC^, 

nNPC^d?/;}) n 

in agreement with the relationship between average linear entropy and delocalization 
of one-particle states found in [27] . 

Remark 3.1. It may be interesting to observe that, for an arbitrary state \tp), it 
is always possible to identify a product basis where 

P,oc(|V')) = l-^E^-J>fc|'l«J-|'- (15) 

k<j 



PLT'\m - .r . .... .u.^^ - (13) 
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To show this, note that the expectation value of an arbitrary traceless normahzed 
ith-qubit observable may be written as {ip\ax^ IV') + ny^ (V'ky 
for some unit vector n*^*' = {nx\ ny\ n^z^ ) . Since the 3-tuple of real numbers ((cri*^), 
{(Jy^), (ci*''}) can be considered a vector in R'^, there clearly exists a direction n*^*-' that 
is parallel to the vector of expectations. One may associate a traceless normalized 
single-qubit operator with the parallel direction for each qubit, a'f' — n*^*-* • 
The mutual eigenstates of the {tri*''} then form a product basis in which Eq. 
holds since the single-qubit operators perpendicular to the {o'i'''} have vanishing 
expectation values. In this basis, the reduced density matrix of each subsystem 
is diagonal. As such, the expression for a state in this basis may be considered a 
standard canonical form which generalizes (non-uniquely) the Schmidt decomposition 
for bipartite systems [551 HHj ■ 

3.2. Generalization to qudit systems 

Eq. (jlSp and, under appropriate conditions. Theorem 3.1, may be generalized to a 
system consisting of n d-dimensional subsystems (qudits). 

To this end, begin by observing that because each reduced density matrix is 
Hermitian, it is always possible (similar to the d = 2 case) to find a product basis 
where each qudit reduced density matrix is diagonal. Let each state in such a basis be 
specified in terms of quantum numbers . . . , where Vi labels a state of the i-th 
qudit and may take any of d possible values. Let \^p) = J^k '^k\k), where fc is a collective 
index ranging over all possible strings of values {vi, . . . ,Vn). The reduced density 
matrix for the i-th qudit may then be expressed as pi = J2v {^k\k =v I'^fcl^) 
and p2 ^ 'E^^{J2k,k'\ki=k',=v,M'^\(^k'\'^)\v^){v,l whereby tr(p2) is the sum over 

all terms |afcp|afc'p such that h = k[. Hence a term |afcp|afc'|^ occurs in the 
sum over different qudits, X]itr(/9f), as many times as Vi = v[. Let fkk' be 
the number of instances in which Vi ^ v[ over all i. This may be considered a 
generalized Hamming distance^ which reduces to the usual one for d = 2. Now 
I]ztr(pf) = J2kk'(^ ~ fkk')\ak\'^\ak'\'^- Using the identity J2kk' = 1, this 

may be rewritten as J^i^^ipf) = " - J2kk' /fcfc' Thus, 

When the value of d is such that a maximal set of (d + 1) mutually unbiased 
bases spanning the state space of each qudit exists, the construction in |30| implies 
the existence of an Hermitian operator basis for the unilocal observables on each qudit 
which is partitioned into (c?+ 1) maximally commuting subsets. Accordingly, the local 
purity may be written as P/oc = X^aii Pa, where Pa is the purity with respect to 
a choice of one maximally commuting set of basis operators for each of the qudits. 
The operators contributing to each Pa uniquely define (up to irrelevant relabeling 
transformations) a product basis, {Ifca)}, where ka is a collective index for the local 
quantum numbers, u^, which label the mutually unbiased bases of each qudit. Since 
for any a and /? \{vi---Va\v}---v^)\'^ = = 1/^" foi" values of {w^} 

and {■y^}, it follows that the product bases {1^^)} are mutually unbiased. 

Recall that Eq. (fTB)) results from summing the squared diagonal matrix elements 
of each reduced density operator in a particular product basis. The stipulation that 
each reduced density matrix is diagonal in this basis ensures that such a sum yields 
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Pioc (after subtracting the trace contribution and proper normalization). If such 
a condition is relaxed, then Eq. (fTB)) states a relationship between the purity with 
respect to an operator basis spanning the diagonal observables of each qudit and the 
components of state vectors along the corresponding product basis. Hence Eq. (jl6p will 
in general hold between each Pa and the state vector components in the corresponding 
basis {Ifca}}. Considerations similar to the ones presented for the qubit case are then 
applicable. Thus, Pioc may be related in general to NPC, and correlations with respect 
to (generalized) Hamming distance in each mutually unbiased product basis. 

Remark 3.2. Interestingly, the expression X]fc<fc' /fcfc' kfePkfc' P be 
interpreted as the expectation value of the (generalized) Hamming distance between 
measurements on two copies of the state \ip). Thus, the local GE, GEjo^ = 1 ^ Pioc, 
may always be written as 

n[d — 1) 

where -F is a Hermitian operator which may be interpreted as the Hamming distance 
between measurements on two copies of the same state in the canonical, state- 
dependent basis in which the subsystem reduced density matrices are diagonal. 

For d such that each qudit may be spanned by each of a maximal set of mutually 
unbiased bases, GEjoc may additionally be expressed in terms of the expectation values 
of Hamming distance between measurements of two copies of the same state in each 
of {d + 1) mutually unbiased product bases, 

2d "^^^ 
GE,„,(|V^)) = — — -(V'|®(^|(^F„) \^)®\^)-d. (18) 

' a=l 



4. Average generalized entanglement of random pure states 

The requirement of Hamming uncorrelation which is responsible for a simple 
relationship between global multipartite entanglement and delocalization is naturally 
satisfied on average by certain classes of random states. 

One such family may be defined, for instance, by taking an arbitrary set of 
normalized probabilities, assigning them at random to basis states in {|fcz)}, and 
giving each component a random phase. The resulting NPCz is determined exactly by 
the set of probabilities, and is the same for all states in the ensemble. The distribution 
of components in the {Ifc^:)} and {jAj,}} bases, and hence the expected value of 
NPCx and NPCj, is determined by the set of probabilities. For particular states of 
the ensemble NPC^; and NPCj^ will fiuctuate around this value. The random assignment 
of probabilities ensures that in the {|fcz)} basis no correlation between component 
products and Hamming distance exists on average. Furthermore, the random phases 
ensure Hamming uncorrelation in the {lA:^:)} and bases also. Thus, ensemble 

averages over many random assignments will yield the relationship in Eq. (|lip . 

In practice, random states generated by uniformly sampling according to the 
invariant Haar measure play an important role, naturally emerging, in particular, 
within statistical descriptions of complex many-body systems such as Random Matrix 
Theory (RMT) [3T]. Results on the expected linear entropy of a subsystem date 
back to early work by Lubkin [32], have been further extended in [33], and more 
recently revisited in the context of obtaining estimates of the expected value and 
variance of the Meyer- Wallach global entanglement J,, and generalizations to other 
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bipartite divisions [31]. Results on the full probability distribution have also been 
established under additional restrictions on the set of states and/or entanglement 
measure [351 136[ 137) . Here, we begin investigating typical GE properties with respect 
to an arbitrary observable set, and show that a simple method allows to calculate the 
expected ft,-purity, P/i, defined in Eq. ([5]). We have the following: 

Theorem 4.1. Let h be any (Hermitian closed) subspace of traceless observables 
on Ti. The expected /i-purity of a pure state sampled uniformly according to the Haar 
measure is given by 

n = E(«-){P,(|^))} = .,^. (19) 

Proof. We first show that the ensemble expectation E{('0|6i|?/')^} is the 
same for any normalized traceless operator spanning h. Let hi — '^Xi\n){n\ 
be a spectral decomposition of bi. Now, E{(V'|6i|V')^} = S '^n^'^IK'^IV')!''} + 
2 ^ A„AmE{| I (m|V') P}. Since by assumption the distribution of \ip) is invariant 

under arbitrary unitary transformations, the expectation E{|{n|'0)|'*} is the same for 
all n, and E{|(n|-0)p|(m|^/')p} is the same for all pairs m ^ n. From the trace and 
normalization conditions, ^ A„ = 0, and -^n = it follows that J2 ^m^n = —N/2. 
Thus, E{{iP\bi\^y} = AfE{|(0|V')|'*} -iVE{|(0|V^)|2|(l|V')n, irrespective of i. 

The value of E{(?/'|&i|V')^} may be determined by using the property that the 
purity relative to the full space of observables equals 1. Since, by Eq. ([2]), KaU — 
1/{N — 1), and {N"^ — 1) linearly independent traceless operators exist, the required 
expectation is 

The expected /i-purity is Ph = Hh X^i which yields the desired result. □ 

Example 1. For a system of n qubits, the local purity of a typical pure state 
averaged over the Haar measure on SU(2") is found to be 

— _ 3n 3 



'N + i N + r 

in agreement with the result for GEioc — Q — {N — 2)/{N + 1) derived in |8]. 

Example 2. As a further application, consider a spin- J system, living in a Hilbert 
space of dimension = 2J + 1, carrying an irreducible representation of SU(2). If 
SU(2) observables are distinguished, the corresponding 5u(2)-purity is 



P.u(2)(|^))-^ E (V'l&dV')', \i J(J+l) -^'^ 

i=x,y,z V ^ ' 

where denote angular momentum operators, and K<iu(2) = ( J + 1)/3J is chosen 
so that Psn(2){\4')) — 1 for angular momentum generalized coherent states [38]. The 
above GE measure may be directly relevant to describe GE generation in a quantum 
kicked top initially prepared in a spin coherent state [12] . In a parameter regime 
corresponding to chaotic dynamics [23j , RMT predicts the long-time asymptotic state 
of the top to be described by a random pure state uniformly drawn according to the 
Haar measure on SU(A^). By the above Theorem, the expected su(2)-purity may then 
be estimated as 

- _ 3 1 
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This coincides with the resuh obtained in [39j by direct integration, and is in exceUent 
agreement with numerical simulations [12] . 

As noticed, for states obeying an appropriate anti-unitary symmetry, the 
components may be chosen real without loss of generality. For random states with 
purely real components, only {N — l)(-j + 1) operators are required to span the space 
of real traceless observables, resulting in 



1 



Kh{N ~-l){N/2 + l) N + 2' 



thereby 



P. = (21) 

where h is now understood as a subspace of purely real observables. 

Example 3. The expected value for the IPR in any given basis for random states 
with purely real components may be found by exploiting the connection between 
IPR and Ph^iag shown in Eq. ([7]). Since (iV — 1) basis operators span hdiag and 
Hh^^^^ = 1/{N + 1), it follows that P/i^.^^ = 2/(A^ + 2). Thus, 

N + 2 ^ ' 

The result given in Theorem 4.1 may also be extended to situations where the 
random states of interest belong to a proper subspace S C H with dim(5) = Ns- In 
general, care should be taken as the basis operators bi need not remain traceless 
and normalized after projection into S. Let 11 be the projector onto S. Then 
nbiU = aib'^ + 13,1, where tr(6';) = 0, and tr(6'- ) = Ng. Now E{(i/'|ai6^ + = 
a^E{(^\b',\^y} + a^PMWKl^)} + P^- But E{{iP\b[\^)} = since E{\{n\^)\^} does 
not depend on n, and b', is traceless. Thus, by Eq. ((20)) one finds 



p^is--^{j^,T.-^ + T.p^)- (23) 

■i i 

Example 4. Consider the average local purity for pure states of the Sz — 
subspace Sq in the state space of n qubits, which have real components when expressed 
in {lA:^)} basis, and are uniformly random with respect the Haar measure on SO(A^o)i 
dim(5o) = Nq. The only singic-qubit observables having non-vanishing expectation 
values for states of this ensemble are cri'-*. Since each o-i'-* is diagonal in the (Ifc^)} 
basis, Hcri'-'n is also diagonal. Furthermore, since every (diagonal) matrix element is 
either -|-1 or —1, it follows that tr((n(Ti*''n)^) — Nq- But because there are as many 
{|fcz)} basis states spanning Sq for which the i-th qubit is as 1, it also follows that 
tr(ncri*'n) = 0. 
Haar measure is 



tr(ncri*''n) ~ 0. Thus, the local purity of a typical real pure state averaged over the 



— 2n 

Ploc\So — 



No + 2 No + 2 

Example 5. A similar method may be followed to obtain the expected 
purity with respect to other subalgebras of qubit observables, in particular algebras 
corresponding to all observables on selected pairs or g-dimensional blocks of spins 
(e.g. bi- local purity P2, tri- local purity P3, and so on). Consider, for instance. 
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the case q = 2, which is relevant to the analysis in [11 . That is, we wish to 
compute P2, over pure states of the Sz — subspace of an n-qubit space, with 
real components in the {|fcz)} basis, which are uniformly random with respect 
the Haar measure on SO(A'o). Since P2 ~ j^J^iPbh, where Pui is the purity 
if the «-th 2-qubit block, it suffices to calculate Pui- The only two-qubit Pauli 
operators which have non-zero expectation values for this ensemble are: tri^'' , ai'^^ , 



(1) (2) (1) (2) 



and a'^^u'^^ 



The trace and trace-norm of the projection of 
each operator into Sq may be found using combinatorial arguments presented in 



yielding: tr(nCTi^''n 



= tr(n4')n) 



0, tr((nai^)n)2) 



tr((nai'^n) 



A, 

(1)„(2)t 



tr((nai^)ai')n)2) 



^0, 



tr(ncri^Vi^^n) = tr(n4^V^^^n) = O, and tT:{{Iial'-> a),^>nf) - ..yy^.uy uy 
((ra-^va) ■ "^^^ coefficients and (3i for the traceless and identity components of 



\{n-2)/2l 

the projection of each operator into Sq may be determined from these values. Thus, 
applying Eq. (^5]) finally yields 



P2\So - g 



No 



A2 


4 / L-2 Y 






^ No \{L-2)/2). 


^ No- 



5. Application to disordered quantum spin chains 

A natural testbed for the above considerations is the study of many-body quantum 
systems. Here, we focus on investigating the relationship between local purity and 
NPC in the eigenstates of a disordered Heisenberg spin chain across a transition from 
quantum integrability to quantum chaos. 

Quantum chaos is generally understood as referring to manifestations of classical 
chaos at the quantum level. Foremost among these is the distribution of energy level 
spacings. As it is by now well established, classically integrable (chaotic) systems 
typically exhibit a Poisson (Wigner-Dyson) level statistics distribution 23J. For 
systems without an obvious classical counterpart, for instance spin chains, the presence 
of a Poisson or Wigner-Dyson level spacing distribution is taken as a phenomenological 
criterion for labelling the system as integrable or, respectively, chaotic. 

In what follows, we shall consider a representative disordered quantum spin 1/2 
system within a class of Heisenberg models in a transverse field which we discuss in 
full generality in [40j . In particular, we choose a one-dimensional quantum spin chain 
described by the following Hamiltonian: 

n J n—1 

H = Ho + H^nt = E J-^^ + I E ^^^^ • '^^'^'^ ' (24) 
1=1 j=i 

where = e + Ssi, e and J are fixed positive numbers. Set are uniform random 
variables within the interval [—d/2, d/2], and open boundary conditions arc imposed. 
Because H commutes with the z-component of the total spin angular momentum 5*^, 
each invariant subspace may be diagonalized independently. We focus on the band 
with no net magnetization, the Sz = subspace. Eigenvalues and eigenvectors have 
been computed numerically for chains of size up to n = 12. This yields N = ( g^) — 924 
as the dimension of the relevant Sz = subspace. 

When J/d — 0, H = Hq is trivially solvable, and for sufficiently small J/d 
perturbation theory is valid. In this regime, the system has Poisson level statistics. 
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When d J, perturbation theory breaks down, and a cross-over from Poisson to 
Wigner-Dyson level statistics occurs. 

Associated with the transition in level statistics, there is a transition from 
eigenvectors which are well localized in the eigenbasis of Hq, to eigenvectors which 
are delocalized and approximately random (Fig. [J). Generally, for fully developed 
chaos, the eigenvectors achieve a distribution of components uniform over the surface 
of an A''-sphere. In systems obeying time-reversal invariance or, more generally as 
mentioned, an appropriate anti-unitary symmetry [23l |41], this is equivalent to a 
Gaussian distribution of eigenstate components in the limit of large N. However, as 
seen in Fig. [1] for this model the states for which NPCz is near the expected value 
for random states of (A^-|-2)/3 [42] have a component distribution which is only 
approximately Gaussian. Furthermore, there is no regime where most eigenvectors 
have an NPCz consistent with the expected value for random states, although this 
value does serve as an approximate upper bound on delocalization. At specific J/d 
values, this model typically exhibits a fairly wide distribution of NPCz. 

Throughout the localized-to-delocalized transition, we examined the relationship 
between NPC2 and local purity for each eigenvector in the Sz — subspace. In 
Fig. El Pioc is plotted against NPCz for each eigenvector using a single random disorder 
realization and four representative values of J/d. At J/d = 0.59, Pioc is averaged over 
each eigenvector between integer values of NPCz and over 100 disorder realizations, 
resulting in a smooth curve which closely fits 

Piocili'))^ 0.032, (25) 

^'^'^ NPCz (It^)) + 12.2 

see inset in Fig. [21 The value J/d = 0.59 is chosen because of the corresponding wide 
distribution of NPCz. Similarly constructed average curves for other values of J/d, 
however, do not show significant differences. 




1 50 200 
NPC 



350 



Figure 1. Eigenvector distribution versus NPCz at J/d = 0.20, 0.59, and 1.0, for 
model Hamiltonian l|24|l with n = 12 spins. Inset: Distribution of eigenvector 
components, a, for eigenvectors with 300 < NPC^ < 316, based on 300 random 
realizations at J/d = 1. The smooth curve is a Gaussian distribution with 
0-2 = 1/924. 
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NPC. 



Figure 2. Pioc versus NPC^ for each eigenvector in the central band using a 
single random realization at four different values of J/d. Inset: Linear fit of Pjoc 
to (NPCz + 12.2)~^ over all eigenvectors and 100 random realizations at J/d = 0.59. 




Figure 3. Dependence of averages on NPC^ at J/d = 0.59. Each is 
averaged between integer values of NPCz and over all eigenvectors and 100 random 
realizations. 



In spite of qualitative agreement, the averaged relationship between Pioc and 
NPCz given in Eq. ([^5]) deviates from the relationship predicted in Eq. under the 
assumption of Hamming uncorrelation between A J — |a^|2|a|P and /^j. This indicates 
that the AJ: do depend non-trivially on Hamming distance in general. In Fig. [3l v4j 
is averaged over each eigenvector between integer values of NPCz and 100 disorder 
realizations at J/d = 0.59. Especially for small NPCz, Aj tends to be larger for smaller 
values of fkj- There is a strong peak in A2 at NPCz « 4, and a less pronounced peak in 
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Al at NPC^ w 15. All other Ay rise gradually. As NPC^ approaches the limiting value 
of {N + 2)/3, the values of the Af become closer to each other, and appear in good 
agreement with the expected value for random states, {N — 3)/[N{N — 1){N + 2)] 
[13] . AIq and Af 2 (data not shown) lie close to A| . 

In the perturbative regime, the dependence of Aj on Hamming distance fkj may 
be understood as a consequence of the two-body form of the interaction. For small 
J/d, the eigenvectors may be expanded in a perturbation series. Starting with an 
arbitrary eigenvector |fc) of Hq, 

\k) ^ Ifc) + V 1?) + C rIO + ••• 

The diagonal contribution az az may be incorporated into the unperturbed 
diagonal energies for the current reasoning. Since the Heisenberg interaction only 
couples eigenstates of Hq which are Hamming distance of 2 away from each other, 
every state {\j)} that appears in the first-order sum is a Hamming distance of 2 away 
from |fc). The states {\l)} that appear in the second-order sum are Hamming distance 
2 or 4 away from Thus, for an eigenstate in the perturbative regime, all of the 
first-order products will contribute to A^. No product larger than second-order will 
contribute to A|, and so on. After the breakdown of perturbation theory, the Aj 
continue to show a dependence on order in perturbation theory for all values of NPCz, 
however, the effect decreases as NPCz approaches the random state value of {N + 2)/3. 

An interesting question is the behavior of the relationship between NPCz and 
Pioc in the thermodynamic limit where n — > oo. Because the relationship between 
local purity and NPCz involves averages over all pairs of basis states of fixed Hamming 
distance, it is reasonable to conjecture that the relationship should become increasingly 
sharp in this limit, provided that the average value of |afcp|ajp for fixed order in 
perturbation theory exists. 

6. Conclusion 

We have quantified the relationship between delocalization as measured by NPC in a 
maximal set of mutually unbiased product bases and global entanglement as measured 
by local purity. In general, the relationship between the two depends on how products 
of state vector components are correlated with respect to Hamming distance - or 
a suitable generalization for higher-dimensional subsystems. Under the condition 
that no such correlation exists, a simple relationship between NPC in each basis and 
local purity is established. For states with certain physically relevant symmetries, 
the number of bases may be reduced. In addition, for each state, there always 
exists a basis in which the local purity is related to NPC in this single basis through 
correlations with respect to Hamming distance. Such analysis yields an expression 
for local entanglement, GEjoc, as the expectation value of Hamming distance between 
measurements of two copies of the same pure state in the state-dependent canonical 
basis where each reduced density matrix is diagonal. 

Distributions of random states under which the assumption of uncorrelation is 
naturally satisfied are also discussed. A simple method to calculate the expected 
relative purity over an ensemble of pure states invariant under the Haar measure is 
introduced, and illustrated in several examples. Lastly, the connection between local 
purity and correlations between products of components is investigated numerically 
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for a disordered Heisenberg spin chain. Because the deviation of the relationship 
between Pi^c and NPCj, from that predicted under the uncorrelation assumption is 
likely a consequence of the two-body nature of the interaction, a similar relationship 
is predicted to hold for any disordered qubit system with two-body interactions which 
has symmetry properties allowing NPC in a single basis to enter the relationship with 
Pioc- For systems without such symmetries, we conjecture that the Pa associated with 
NPC in the eigenbasis of Hq will still provide the main contribution to Pioc, until the 
eigenvectors are maximally random. Thus, the relationship between Pioc and NPC in 
the eigenbasis of Hq may be generic to all disordered many-body systems. 

As a general remark, we also expect the correlation between products of 
components and Hamming distance to be relevant to other entanglement measures. 
For instance, the n-tangle is written as a sum of products of pairs that are Hamming 
distance n-apart. Thus, we conjecture that this characteristic structure may be 
important for the study of entanglement properties across a localized-to-delocalized 
transition and across quantum criticality in many-body systems. 
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